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Abstract 

We demonstrate through exact solutions that a spin bath leads to stronger (faster) dephasing of 
a qubit than a bosonic bath with identical bath-coupling spectrum. This difference is due to the 
spin-bath "dressing" by the coupling. Consequently, the quantum statistics of the bath strongly 
affects the pulse sequences required to dynamically decouple the qubit from its bath. 

PACS numbers: 03.65.Yz, 03.67.-a, 03.65.Ud, 03.65.-w 



*Electronic address: dasari@weizmann.ac.il 



1 



I. INTRODUCTION 



An essential requirement for realizing fault-tolerant quantum computation is the ability 
;o protect the information encoded in the qubits from leaking into the environment (bath) 
1|. The stronger the qubit-bath coupling, the faster is the loss of information (decoherence) 
2] . The desirable but often unachievable goal is to completely decouple the system from the 
bath. Yet, even if they cannot be completely eliminated the adverse effects of decoherence 
suppressed using control pulses at a rate faster than the inverse memory (correlation) time of 
the bath j^, 4 ] . For a common kind of decoherence known as pure dephasing [3| the problem 
can be solved exactly in certain cases jsj. If the bath is taken to be bosonic, i.e., composed 
of harmonic oscillators, but otherwise general, its control can be analyzed by a second- 
order non-Markovian master equation [4]. This master equation has yielded a formula for 
optimal dephasing control [6] which has been extended to arbitrary forms of decoherence 
Interestingly, similar results are reproduced when the source of dephasing is assumed to be 
classical noise, whose power spectrum coincides with the coupling spectrum of the quantum 



bath 



There has been a growing interest in decoherence control from fermionic (spin) baths 



llj as they form the main source of decoherence for solid-state qubits 12J. Decoherence 



in these systems has generally been treated by classical models: for example, in the quasi- 
static approximation I13I. Il4|. only the (effective) magnetic field distributions caused by the 
bath spins determine the dephasing rates and the optimal spin-echo times. 

In general, dynamical control methods set out to affect the qubit-bath dynamics on a 
;ime-scale much shorter than the bath correlation (memory) times known as the Zeno re gim e 



15] . On this time scale, the qubit-bath interaction dynamics is coherent and reversible 15]. 



The control then moves the qubit's frequency away from the harmful bath modes, thus 
approaching dynamical decoupling |3|. When the control is instead applied in the anti-Zeno 
regime, i.e., on the time scale comparable to the correlation time, it can drastically enhance 
the decoherence rate [15J, H6( . Hence, the knowledge of the cross-over from the Zeno to the 
anti-Zeno regime is crucial for designing pulse sequences to control decoherence. 

The foregoing studies leave us with some important questions: Is the knowledge of the 
bath (or noise) spectrum sufficient to design optimal control ? Does one need a rigorous 
quantum analysis of the system-bath dynamics or can it be substituted by classical-noise 
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analysis {5,6]? How does the bosonic or fermionic nature of the bath constituents influence 
the time scales of dephasing and its control ? Recent works have addressed the differences 
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Here 



between the spin and bosonic baths using second order master equations 
we exactly analyze pure dephasing and its control, exactly for spin and bosonic baths and 
contrast, the two, showing the limitations of second-order master equations for each bath. 

We point out that pure dephasing of a qubit coupled to a spin bath, cannot in general be 
treated within the classical noise description (e.g., by assuming an effective magnetic- field 
distribution 

HQ) 

as opposed to dephasing caused by a bosonic bath. We demonstrate 
that even if the coupling spectra of the two kinds of baths are taken to be identical, the 
resulting dephasing is drastically different. In particular, we show that the Zeno to anti- 



Zeno crossover is not solely determined by the bath coupling spectrum [4|, [15| , but strongly 
depends on the quantum statistics of the bath constituents. Correspondingly, the same 
pulse rate is shown to fall within the Zeno (decoherence-suppressing) regime for one bath 
and within the anti-Zeno (decoherence-enhancing) regime for the other. 



II. QUBIT DEPHASING 

The coupling of the qubit (S) to the spin-bath (cr) is described by the qubit-spin (spin- 
spin) Hamiltonian 

H Sa = A(t)S* +5> fe <x* + + a k ). (1) 

k k 

The first term represents the control Hamiltonian for the qubit (S) with strength A(t), 
Uk represents the (Zeeman) energy of the k th bath spin (ak), and r\k is the corresponding 
coupling strength of the spin to the qubit. The control and the bath act on perpendicular 
Bloch-sphere axes of the qubit x and z respectively. 

Similarly, when the qubit is coupled to a bosonic bath with identical bath coupling 
spectrum, its evolution is governed by the spin-boson Hamiltonian ^ 

H Sb = A(t)S x + J2^blh + J2^kS z (bl + b k ), (2) 
k k 

where b, W represent the creation and annihilation operators of the k th boson (oscillator). 
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A. Pure dephasing without control 



(i) Spin bath: In the absence of control, the system-bath evolution leads to pure dephasing 
of the qubit. In contrast to dephasing by a bosonic bath (App. A) the Magnus expansion 
in the case of a spin bath does not truncate to any order (App. A). Hence, obtaining the 
time evolution operator by this procedure is quite cumbersome. 

In the S z basis of the system (|+), |— )), the Hamiltonian of a qubit in a spin bath Eq. 
(GQ) with A(t) = 0, can be rewritten as 

H Sa = h+\+)(+\ + h;\-)(-\, 

Hf = "A ± E + a ^ ( 3 ) 

k k 

where are the bath operators. The simple form of makes it easy to diagonalize, so 
that we obtain a closed-form equation for the time-evolution operator 

uAt)=[u+(t)\+)(+\ + u;(t)\-}(-\] (4) 

where 

k 

Here the renormalized eigenfrequencies of the bath spins are 5 k = ^oo\ + rjf,. As shown 
below, this renormalization would lead to faster dynamics for a spin-bath than a bosonic- 
bath. 

As the above dynamics leads only to pure dephasing of the qubit, the populations of the 
|+) and |— ) states are unchanged, while the off-diagonal elements Pg(t) dynamically evolve 
in time as 

pt(t) = e- tr °Vpt(0), IV(f) = -~ln {Tr B [U+(t)pAW B (t)}) . (6) 

where p CT (0) is the initial state of the bath spins. 

We consider an uncorrelated initial state of the bath spins p a = exp(— BHb)- In the 
energy basis, the state of a k th bath-spin can be written as p a ^ = Pk\ e )( e \ + (1 ~ Pk)\g)(g\, 
where pk = exp(— /3o;/ c )/2 cosh/Jo^. Since the trace of bath operators obeys the simple 
symmetry rule Tr[U^ k (t)\e) (e\U~ k (t)} = Tr[U£ k (t)\g) (g\U~ k (t)], one immediately finds that 
the dephasing rate is temperature independent and is given by 

rv(i) = ^ E ln ( x - 2 f sin2 6kt ) ■ ( ? ) 
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COS Sutl + 



i sin Skt 



(u k a z k ± n k a x k ) 



(5) 



(ii) Bosonic bath: In the spin-boson model one obtains a closed-form expression for the 
time-evolution operator using the Magnus expansion, which truncates to second order as 
the commutator of the bosonic creations and annihilation operators is a c— number (App. 
A). Equations (jlJ-(JB]) hold true even for a bosonic bath, if are replaced with appropriate 
displacement operators for the oscillators, 

U?(t) = J] U Ut) = e«) 6 *-^\ (8) 

k 

where a = (l—e tuJkt )/uk- Here the eigenfrequencies of the bath particles are not renormalized 
(dressed) by the coupling unlike for the spin-bath. 

Substituting the expression above in Eq. (jSJ), we find that a bosonic bath would lead to 
temperature-dependent dephasing rate 

1 n 2 

T b (t) = - V 4 coth (3u k sin 2 u k t. (9) 
1 k u k 

At low temperatures one always finds that 

r a {t) > r b {t) (io) 

i.e., a spin bath leads to more harmful effects than a bosonic bath. In Fig. 1(a) we compared 
these rates for a Lorentzian bath. 

It can be seen from Eq. (j7]) that the spin-bath leads to complete dephasing at certain 
times if there exists even a single mode whose coupling strength satisfies % > Uk- This 
is quite dramatic as compared to the oscillator bath which cannot lead to such complete 
dephasing even if all the bath modes satisfy the condition rjk > ujk- 

The other significant difference is the faster dynamics that the spin-bath induces due to 
its renormalized frequencies 



S k = yu 2 k + Vk>Uk- (11) 
Because of this faster dynamics, despite a similar bath coupling spectrum, a given time 
scale in the Zeno regime for spin-boson dynamics can lie in the anti-Zeno regime for the 
spin-spin dynamics. This can make a huge difference to the control schemes which strongly 
rely on the appropriate time scales for operation. In addition, the effect of renormalized 
frequencies becomes more evident in situations where the dominant coupling of the qubit is 



to the zero-frequency modes (1// noise [19j). In such cases the bosonic bath leads to slow 
exponential decay even at short times (Markovian dynamics), while the spin-bath can still 
lead to faster Gaussian decay (non- Markovian dynamics). 
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FIG. 1: Coherence dependence on quantum-statistics for qubit coupled to a Lorentzian 
bath. (a). The free evolution of the qubit coherence pt~(t) is shown to be faster when the bath 
consists of are either spins (solid-blue line) than when it consists of oscillators (dashed-green line). 
Time is in coupling strength units ry _1 . In the inset shown is the evolution for longer times, (b). 
We have compared the relative change (enhancement) of qubit coherence in the presence of control 
with respect to its free evolution for spin and oscillator baths, under periodic 7r-pulses applied at 
a rate r = T/n, where T is the total time and n is the total number of pulses. 

B. Dephasing control for Spin and Boson baths 



We consider brief control pulses along x at intervals r as a means of dynamically de- 
coupling the qubit from the bath. The simplest control that can significantly reduce the 
dephasing rate can be realized by applying frequent brief 7r-pulses in a Bloch-sphere direction 
(x) orthogonal to the bath- interaction (z). 

(i) Bosonic bath: A closed form equation for r#(£) can be obtained for periodic 7r-pulses, 



t k u k 



sm 2 (f) sin 2 '^±^ 



(12) 



sin 2 (^) sin 2 (*f) 

For bosonic baths one can obtain exact analytical expressions for the dephasing rates even 
for aperiodic 7r-pulses. In particular with high-frequency cutoff, an aperiodic pulse interval 
r n = Tsin 2 (jrf), where T is the total time and N is the total number of pulses, was shown 
to be the optimal pulse sequence for decoherence control [5|. Yet, this pulse sequence may 



not be optimal if there is a constraint on the pulse number (energy constraint) 3|. Since 
the minimal pulse interval is determined by the total number of pulses N, the choice of N 
determines whether these pulses give rise to the Zeno or anti-Zeno effects {4]. 
(ii) Spin bath: The form of T a (t) in Eq. ([7j) makes it difficult to obtain a closed-form 



expression when brief 7r pulses are applied at aperiodic intervals r n . Yet, the nature of 
spin- 1/2 nature of the bath constituents allows one to exactly diagonalize the system-bath 
dynamics numerically for any number of bath particles in the presence of control. The 
control dynamics is determined by multiplying many 2x2 matrices, which is evident from 
the structure of the unitary evolution operator: 

Ua(t)=[---U+(r 3 )U-(T 2 )U+(n)] \+)(+\ + [---U-(r 3 )U+(r 2 )U;(r 1 )] |-)(-| (13) 

A simple closed-form expression for T a (t) when 7r-pulses are applied periodically at in- 
tervals r, (see App. B) is found to be 

= ^ E ln f 1 " 2^ fe (r)| sin 2 S k t\ , F^(r) = 
k k / 

where 

fc = cos _1 (l - —± sin 2 5 k t). (15) 

As discussed earlier, the faster dynamics induced by the spin-bath modifies the crossover 
from the Zeno to anti-Zeno regime for a given bath spectrum. We have shown this explicitly 
in Fig. 1(b) and Fig. 2 for the case of a qubit coupled to a Lorentzian bath spectrum, 
where the periodic 7r-pulse control reduces the dephasing rate for a bosonic bath, while the 
same control enhances the dephasing rate for a spin bath. On the other hand, as the control 
pulse interval becomes shorter, the qubit dephasing is more strongly suppressed in the case 
of a spin-bath than in the case of an oscillator bath (see Fig. 2). This contrasting behavior 
stems from the faster dynamics in the presence of a spin-bath, due to dressing of the bath 
frequencies by the coupling strengths. 

III. WEAK COUPLING LIMIT 

n 

For boson baths the second-order (time-local) non-Markovian master equations [4J are 
known to be adequate for evaluating the qubit 's decoherence rates and its modifications 
in the presence of control. This is true even in the limit of strong system-bath couplings, 
due to the commutator of the bosonic creation and annihilation operators being a c-number 
(App-I). As the spin-bath operators do not allow for this simplification it is instructive to 
discuss the limitations of the second-order analysis in the system-bath coupling strength. 



sin n<pk 



2 sin 0^ 



(14) 
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FIG. 2: Quantum Zeno (QZE) to anti-Zeno (AZE) cross over for qubit coupled to 
a Lorentzian bath under periodic 7r-pulses. The relative change (enhancement) of qubit 
coherence in the presence of control with respect to its free evolution for spin (left) and oscillator 
(right) baths, is plotted as a function of the pulse interval r (in units of inverse bath correlation time, 
r) and pulse number n. The transition from QZE (coherence enhancement) to AZE (coherence 
suppression) for bosonic baths is exactly at the bath correlation time 1/T, while for spin baths it 
is well within the bath correlation time. 

(i) Spin bath: The differential master equation (ME) for the density operator of the system 
Ps(t), to second order in the system-bath coupling, is given by [l] 

ps(t) = -i[H s (t),ps(t)) + f dt'mt-t')[S(t)p s (t),S(t')}+h.c.}. (16) 

J o 

The bath response function $(£ — £') = ^fe^fc [cos(u>k(t — t')) + i(ofc) sva{u^{t — t'))] and 
S(t) = U s (t)S z Ul(t), where U s (t) = exp( /„* A(t')dt'S x ). In the absence of control (A(t) = 
0), the second order analysis leads to pure dephasing, at a rate 

T^ IE {t) = Y j %{uj k t-smLo k t). (17) 

On comparing this expression with the exact form for a spin bath, Eq. one can see that 
the only similarity is the temperature-independent dephasing and 0(t 2 ) dynamics. Clearly, 
the second-order analysis overestimates the decay rate and does not capture the effect of 
renormalized system-bath couplings, as obtained from the exact analysis. 
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Under periodic 7r-flip control, Eq. (|T4|) is replaced by 

k ^ 

where the expression for F^ k is given in Eq. ( IT2|) . 

The validity condition of second-order analysis for spin baths is stringent, i.e., 

Vk/wk < 1 (19) 

(ii) Bosonic bath: For a bosonic bath, the second-order analysis leads to an identical de- 
phasing rate for both strong and weak coupling, as that obtained from the exact analysis, 
viz., 




Similarly, under periodic pulse pulse control, the expression for r#(£) is identical to that 
obtained under exact analysis (see Eq. ( TT2|) ). 

For the second-order analysis to be valid, coupling to a bosonic bath only needs to satisfy 

r] k /uj k < 1. (21) 

IV. CONCLUSIONS 

In conclusion we have shown that the temperature-independence of the dephasing rate 
and the "dressing" of the bath frequencies by their coupling strengths render the dependence 
of the dephasing of a qubit by a spin bath very different from that caused by the well studied 
bosonic (oscillator) bath. Consequently, the quantum statistics of the bath determine the 
crossover from the Zeno to the anti-Zeno regime. This in turn affects the decoherence control: 
pulse timings required for Zeno-like dynamical decoupling depend on whether the bath is 
composed of spins or oscillators. Since pulse errors commonly increase monotonically with 
the number of pulses, there is a lower bound on r in any experimental scenario. The useful 
pulse intervals for coherence control will therefore crucially depend on whether the qubit is 
coupled to a spin or bosonic bath. In particular, short-interval (high-rate) pulses can cause 
the leakage of the population out of the qubit subspace to higher levels of the system Q]. 
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Appendix A: Magnus Expansion 



U(t) = T<_exp jT dt'H^t') 



To obtain a closed-form expression for the time-ordered unitary operator 

(Al) 

we resort to the Magnus expansion [20|] of the exponent of U(t) = exp(p,(t)). The first few 

terms of the expansion are 

ft -| ft r t\ 

n(t) 



f H I {t 1 )dt 1 -\ f [ [H^Hifo^dtrfh 
Jo * Jo Jo 



+IJ J 1 J* dhdhdt^Hrih), [Hj(t 2 ), Hj(t 3 )] + [[Hjih), H z {h)], H^t,)}) +(-A2) 

(i) Bosonic baths: Let us consider Eq. (1) of the main text, and set A(i) = 0. In the 
interaction picture this yields 



H T (t) = e iHBt H T e- iHBt = S z J2vk(e~ iu>kt b k + e~ iuikt b\). 



(A3) 



For bosonic bath operators, the commutator of the interaction Hamiltonian at two different 
times is a C-number function in the bath operators: 



[Hi{t), #,(*')] = -2z Vt sinw fc (t - 0- 



The fact that this commutator is a C-number implies that only the first two terms of the 
expansion are non-zero. Now, the closed-form equation for the time-evolution operator takes 
the simple form 



uit) = 5> ± (*)i±x±i 



U±(t) = exp 



The coupling to the bath determines the functions 



tf{t) ± £ (a k (t)b{ - ai(t)h 



f(t) = -^vK^kt - smu) k t)/u)l, a k (t) = r] k 



1 - e 



(A4) 



(A5) 
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(ii) Spin baths: Similar to the above analysis, in the interaction picture we get 

Hi{t) = e iHst H ie - iHBt = S Z J2 Vk{e- iUkt ^ + e"^V+). (A6) 

k 

For spin bath operators, the commutator of the interaction Hamiltonian at two different 
times is an operator: 

k 

The fact that this commutator is a not a C-number implies that all terms of the expansion 
are non-zero and it is quite complex to obtain a closed-form equation for the time-evolution 
operator using magnus expansion. 

Appendix B: Periodic 7r-pulse control for spin baths 

For 2n periodic 7r-pulses, Eq. f[T5j) gets simplified to 

17(f) = [Uz(T)U+(T)] n \+)(+\+ [UUr)Un(r)] n \-)(-\. (Bl) 

Hence a closed form expression for U(t) can be obtained by diagonalizing the matrix 
[UZ(r)U+{T)] n . Since 

[U B (r)UUr)] n = J] KbKbT , (B2) 

n, 

it is sufficient to diagonalize the effect of one bath spin on the qubit after 2n pulses. 

Since the eigenvectors of a unitary matrix B U^ B are orthogonal, evaluating arbitrary 
powers of this matrix becomes quite simple. 

tt- tt+ ( 2r u l ~ Vk ■ 2 c A t 2ico k sm5 k t . v 

U k,B U k,B = I cos °kt p — sm °kt J H p \°k cos °ktcr k - Vk sm b k ta%) 

= KH)(vt\ + KK)(%\- (B3) 

The corresponding eigenvalues and eigenvectors are given by 



x k = x k ±iJl-x 2 k , 



H) = -J=^[\0)k T ia k \l) k ), a = y * . (B4) 

V 1 + a i V 1 - X i + Vk 



where 



x k = 1 - -TT sin2 $kt, Vk = — sin25 fe f. (B5) 
K Vk 
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Now using the fact that {v£\v k ) = and (v k \v k ) = 1, 

[u^bKbT = (K) n K)(vt\ + (K) n K)(v k 



(B6) 

Substituting this in Eq's flB2[) . (IB 1 1) and performing the trace over bath degrees of free- 
dom, one can obtain the dynamically modified dephasing rate of the qubit, to be 

smnd>k 

( myL-M- k [T)-^sm-o k ij, * a ] k {T) = 

where 



Fs(t) = ^ E ln ( X - 2^ fc (r)| sin 2 8 k ?j , F™ k 



2sin< 



(B7) 



cos 



2uj 



-sin 2 5 k t). 



(B8) 
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